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We present a model for molecular materials made up of polar and polarizable molecular units. 
A simple two state model is adopted for each molecular site and only classical intermolecular in- 
teractions are accounted for, neglecting any intermolecular overlap. The complex and interesting 
physics driven by interactions among polar and polarizable molecules becomes fairly transparent in 
the adopted model. Collective effects are recognized in the large variation of the molecular polarity 
with supramolecular interactions, and cooperative behavior shows up with the appearance, in at- 
tractive lattices, of discontinuous charge crossovers. The mf approximation proves fairly accurate in 
the description of the gs properties of MM, including static linear and non-linear optical suscepti- 
bilities, apart from the region in the close proximity of the discontinuous charge crossover. Sizeable 
deviations from the excitonic description are recognized both in the excitation spectrum and in 
linear and non-linear optical responses. New and interesting phenomena are recognized near the 
discontinuous charge crossover for non-centrosymmetric clusters, where the primary photoexcitation 
event corresponds to a multielectron transfer. 

PACS numbers: 78.20.Bh, 78.40.Me, 71.35.-y, 42.65.An 



I. INTRODUCTION 

The promise of molecular materials (MM) for advanced 
applications is impressive: not only old devices can be 
fabricated by light-weighted, flexible and cheap materi- 
als, but a new generation of devices is likely to appear 
to drive a molecular electronic and/or photonic revolu- 
tion. Large non-linear responses to external perturba- 
tions (including electrical and optical fields) are required 
for applications and functional MM are often based 
on 7r-conjugated molecules and polymers. Structure- 
properties relationships are well understood at the molec- 
ular level at least for a few families of molecules and 
the guided synthesis of molecules with specified behavior 
has already been demonstratedi^ But as the synthetic 
ability evolves from the molecular to the supramolecu- 
lar level^ interpretative tools are required to guide the 
synthesis at all stages, and supramolecular structure- 
properties relationships must be devised. 

Extending our comprehension of the properties of MM 
from the molecular to the supramolecular level is a chal- 
lenging task since non-additive, collective behavior ap- 
pears in MM as a results of intermolecular interactions. 
Particularly important collective effects are expected in 
functional MM based on large (and largely polarizable) 
7r-conjugated molecules, as it was suggested by Mc- 
Connell, 40 years ago,- and as it was more recently under- 
lined with specific reference to MM for non-linear optical 
(NLO) applications i 

The mean-field (mf) approach is a simple approxima- 
tion scheme to describe the ground state (gs) of MM. 
Within mf the problem of N interacting molecules fac- 
torizes into N effective molecular problems and collectiv- 
ity is introduced by the local self-consistent fields gener- 
ated by the surrounding at each molecular site. Recently 
the mf approximation was implemented into a quantum 



chemical calculation of the polarization of crystals and 
films of large conjugated molecules&2iiS The charge dis- 
tribution on polarizable molecular sites is strongly af- 
fected by local electric fields, and prominent collective 
behavior is recognized in the very large on-site charge 
redistribution in crystals of quadrupolar molecules So 
far the approach was not applied to crystals or films of 
polar molecules, where even larger effects are expected. 

The mf approximation does not apply to excited states, 
that in MM are usually described within the excitonic 
approach, a very popular approximation scheme origi- 
nally developed for weakly interacting MMi^i^ In this 
limit sizeable collective phenomena are expected only if 
intermolecular interactions mix up degenerate (or quasi- 
degenerate) states. Within the excitonic model, then, 
any interaction among states with a different number of 
excitations is disregarded. The mf approximation sup- 
ports collective behavior, but suppresses any correlation 
of the electronic motion on different molecules. The exci- 
tonic approximation instead describes collective and cor- 
related excited states^ But, not accounting for the mix- 
ing of states with a different number of excitations, the 
excitonic model disregards the molecular polarizability, a 
serious limitation for MM of interest for advanced appli- 
cations. Moreover, the excitonic model relies on the def- 
inition of molecular (local) ground and excited states, as 
states relevant to the molecule embedded in the material, 
but no hint is given on how to derive these effective states 
from the molecular and supramolecular structure^iS^ 

In clusters of non-polar molecules leading supramolec- 
ular interactions are represented by electrostatic inter- 
actions between transition dipole moments: the result- 
ing J-term in the excitonic model describes exciton 
hopping between different sites^ The same interaction 
is also responsible for the so-called non-Heitler-London 
term: 12 this term, neglected in the excitonic approxima- 
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tion, mixes up states whose exciton number differs by 
two units. A recent detailed study has demonstrated that 
this term has indeed only minor effects in aggregates of 
non-polar chromophores 

The simple excitonic model derived for materials made 
up of non-polar molecules is often adopted, with no 
further scrutiny to describe materials based on polar 
molecules. In this case, however, additional interactions 
appear due to the finite permanent dipole moments of the 
molecular units in both the ground and excited states. 
In particular, the interaction between permanent dipole 
moments on different molecular units introduces in the 
Hamiltonian an exciton-exciton interaction term. This 
term conserves the exciton number and therefore survives 
in the excitonic approximation. The role of the exciton- 
exciton interaction in promoting bound-exciton states 
(and specifically biexcitons) in layers of polar molecules, 
was investigated in ref. 1 151 The formation of exciton 
strings was also discussed in ref. ^3 as a consequence of 
exciton-exciton interactions in models for linear aggre- 
gates. This same concept was applied in a different con- 
text, to explain anomalous spectral features in differen- 
tial transmission spectra of organic charge-transfer (CT) 
salts with a mixed donor-acceptor (D-A) stack motifji&iS 
leading to the experimental confirmation of bound exci- 
ton states. 

The excitonic model was developed to describe MM 
in the limit of weak interactions, i.e. for intermolec- 
ular energies much smaller than intramolecular excita- 
tion energies>ii*i£ Electrostatic interactions among polar 
molecules are instead large, and can easily be of similar 
magnitude or even larger than the small excitation en- 
ergies typical of largely polarizable molecules. The fail- 
ure of the excitonic model to describe interacting polar 
and polarizable molecules is therefore hardly surprising. 
However, since current understanding of optical proper- 
ties of MM generally relies on the excitonic picture, it is 
important to explicitly discuss the limits of this approx- 
imation. We will show that the properties of MM based 
on polar and polarizable molecules are quantitatively as 
well as qualitatively different from the sum of the molec- 
ular properties. New and interesting physics appears in 
these materials due to supramolecular interactions, and 
it cannot be appreciated within standard approximation 
schemes. 

In this paper we discuss the gs and the excitation spec- 
trum of a simple model for clusters of polar and polar- 
izable chromophores. Each molecular unit is described 
in terms of the Mulliken models This two state model 
was originally proposed to describe D-A CT complexes in 
solution^ Later it was adopted to model NLO responses 
of push-pull chromophores j2i an interesting class of tt- 
conjugated molecules with electron D and A end groups. 
In this context the Mulliken model, extended to account 
for the interaction with internal vibrations and/or with 
solvation degrees of freedom, offered an accurate descrip- 
tion of low-energy spectral properties of solvated push- 
pull chromophores Intermolecular interactions are 



introduced in terms of classical electrostatic forces, any 
intermolecular charge resonance being disregarded (zero 
intermolecular overlap approximation) £ The resulting 
model quite naturally applies to MM based on push- 
pull chromophores. By choosing different geometries for 
the cluster, the model describes different kinds of materi- 
als, like, e.g., aggregates, Langmuir-Blodgett (LB) films, 
functionalized polymers, and molecular crystals. 

In a different context, one-dimensional (ID) models 
for aligned polar and polarizable units, only interact- 
ing via electrostatic forces offered a first description of 
the neutral-ionic (N-I) phase transition observed in CT 
crystals with a mixed stack motif>2i This model is fairly 
crude, since it reduces a stack of overlapping molecules 
(...D-A-D-A....) to a collection of non-overlapping dimers 
(...D-A D-A....). In spite of that, a mf description of 
the relevant gs offered a first picture for the N-I transi- 
tion, and specifically for the crossover from a continuous 
N-I interface to a discontinuous transition at large elec- 
trostatic interactions*^ More recently the same model, 
within the excitonic approximation, yielded an interest- 
ing picture for the excitation spectrum of mixed stack 
CT salts with a largely N gsi2J£ 

II. THE MODEL 

The Mulliken models was originally proposed to de- 
scribe charge resonance in isolated D-A complexes on 
the basis of the two limiting neutral (N) and ionic (I) 
states: \DA) and \D + A~), respectively. The same model 
describes intramolecular charge resonance in push-pull 
chromophores where the same two basis states describe 
the limiting apolar and charge-separated (zwitterionic) 
structure^ as schematically shown in Fig. 1 for a typical 
chromophore. Following Mulliken, we disregard all ma- 
trix elements of the dipole moment operator, except the 
dipole moment of the charge separated, ionic (I) state, 
(D + A~ \jl\D + A~) = ea, where a measures the distance 
between D and A molecules in CT complexes or an effec- 
tive molecular length in push-pull chromophores. If 2zo is 
the energy gap between the two basis states, and — v2i 
is the mixing matrix element (in the following we will 
fix y/2t as the energy unit), the Mulliken Hamiltonian is 
written as: 

H = 2z a p - V2ta x (1) 

where p — a z + 1/2 is the ionicity operator, measuring 
the degree of CT from D to A, and a x / z are the Pauli 
matrices. 

The solution of this Hamiltonian is trivial. The result- 
ing ground and excited state are linear combinations of 
the two basis states, and the amount of mixing is defined 
in terms of a single parameter: 

"-"VlTT (2) 
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therefore model each molecule in the I (\D + A~)) state 
as a segment of length a bearing +/— charges at the 
D/A ends, and then calculate Vij based on unscreened 
Coulomb interactions. Specifically, for the three ID lat- 
tices sketched in Fig. ^we have: 



A P- 



D p+ 



Vij = {±iy- j 2vw 



(5) 



where the +/— sign refers to A/C lattices, and 
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Vi4 = 2vw 
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FIG. 1: Upper panel: schematic representation of a push-pull 
chromophores in terms of the two limiting neutral (N) and 
charge-separated, ionic (I) forms. The average charge sepa- 
ration in the ground state, p, is proportional to the molecu- 
lar dipole moment, shown in the figure as an arrow. Lower 
panel: the three one-dimensional lattices considered in this 
paper. Each molecule is represented by an arrow, showing its 
ground-state dipole moment. 



as follows: 



for B lattices. In all cases dy is the distance between sites 
i and j, in units of the distance r between adjacent sites, 
v = e 2 /a is the magnitude of the electrostatic interaction 
between two charges at distance a, and w = a/r. In the 
following we consider v =1 or 2, that, for typical y/2t — 
1 eV valuesj22i2ii correspond to molecular lengths a ~ 14 
or 7 A, respectively. Results will be presented for w in- 
creasing from (the gas phase limit) up to a few units 
(intermolccular distance corresponding to a fraction of 
the molecular length). For B-lattices w — 1 represents 
an upper bound. For clusters of N molecules the above 
Hamiltonian is written on the 2 N basis obtained from 
the direct product of the two basis states located on each 
molecular site. The eigenstates are obtained by exact di- 
agonalization for systems with up to 16 site, by exploiting 
translational symmetry. 



\G) = ^T~p\DA)+^-p\D+A-) 

\E) = Jp\DA) - ^Y~ P \D+A-) (3) 

p is the gs expectation value of the charge operator 
and is proportional to the gs molecular dipole moment: 
(G\p\G) = eap. 

We describe a MM made up of polar and polarizablc 
molecules as a collection of Mulliken-molecules (i.e. of 
molecules described by the Hamiltonian in Eq. Q) in- 
teracting via classical electrostatic forces. The relevant 
Hamiltonian reads: 

H = Y^ ( 2z o/5 t - V2ta Xt i) + ^2 ViiPiPj ( 4 ) 

i i>j 

where i,j count the molecular sites, and the index on p 
and a x specifies them as operators working on the corre- 
sponding site. Vij measures the electrostatic interaction 
between two fully-I molecules on sites i and j. Different 
choices are possible for Vij: dipolar interactions are un- 
realistic for typical push-pull chromophores whose length 
(~ 10 A) is similar to the intermolecular distances. We 



III. THE MEAN FIELD AND EXCITONIC 
APPROACHES 

In order to better appreciate the physics of interacting 
polar and polarizable molecules, we will compare exact 
solutions of the Hamiltonian in Eq. with those ob- 
tained within the mf and excitonic approximations. Both 
approximations are based on the assumption that in the 
gs the electronic motion on different sites is uncorrelated, 
so that the gs describes a collection of molecules each one 
in the local (molecular) gs: |0) = \G\G2....Gi....), where 
|0) represents the mf gs wavefunction for the cluster, and 
\Gi) is the gs wavefunction for the i-th site. For a cluster 
of non-overlapping Mulliken molecules the local gs, |Gj), 
necessarily retains the same form as in Eq. Q , but with p 
depending on the surrounding. In any case, for any given 
p, \Gi) and \Ei) define a couple of local mutually orthog- 
onal wavefunctions, and the 2 N basis functions obtained 
from the direct product of these local states is a complete 
basis for the cluster. We define the Pauli matrix opera- 
tors working on the local \d) and \Ei) basis, S z / Xi i, as 
linear combinations of the Pauli operators, <7 Z / X .i, work- 
ing on the original \DA) and \D + A~) basis, as follows: 
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S z ,i = (1 - 2p)a Zti + 2y/p(l - p)<J Xt i 

(1 - 2p)a x , 



Hamiltonian for our model by fixing p as to impose the 
exact vanishing of this term: 



S x ,i = — 2-\/p(i - p)a 



(7) 



To make contact with standard notation^ we construct 
the following hard-core boson operators: 



s x ,i = (ty + h) 

S zA = l-2b% 



(8) 



These operators have a very clear physical meaning: 
working on the vacuum state |0), defined as the state 
where all sites are in the local \Gi) state, b\ switches the 
i-th site to the local \Ei) state, and hi = SjSj returns the 
and 1 value for the |G,) and |JSj) state, respectively. 
In terms of these operators the Hamiltonian in Eq. Q 
reads: 



H = H m f + H ex + H u 



(9) 



where: 



hi (10) 



H, 



H ™f = [ 2 ( X ~ 2 -°)( z o + Mp) + 4v/p(l - p) 

i 

+ £ [2y/p(l- P )(zo + Mp) - (1 - 2p)] (SJ + h) 

i 

]T Vy [p(l ~ p)(&& + Ml) + (1 - 2pfh i h 3 ]ll) 

Y,Vij [p(i-p)(6iSi+sts|) 

+2(l-2p)Vp(l-p)(fo! + ^] (12) 

and M = E l>3 V^/JV. 

Three terms in the Hamiltonian in Eq.@ have been 
conveniently grouped: the first one, H m f, collects on-site 
terms, and, as it will be discussed below, sets the basis 
for the mf description of the gs. Both H ex and H uex 
work on two sites: H ex conserves the total number of 
excitations and hence enters the excitonic Hamiltonian, 
whereas H uex mixes up states differing by one or two 
excitation, it represents an ultraexcitonic term that will 
be disregarded in the excitonic approximation. 

The transformation in Eq. (7J| holds true for any < 
p < 1, and the total Hamiltonian in Eq. © is not affected 
by the specific choice of p. Of course, its partitioning 
into the three terms instead varies with p, and different 
excitonic approximation can be derived from the same 
model Hamiltonian, depending on the specific choice of 
p, or, equivalently, on the specific choice of the local \d) 
and \Ei) wavefunctions. It is important to realize, how- 
ever, that the second term in H m f does not conserve 
the exciton number and must be killed in the excitonic 
approximation. We therefore define the best excitonic 



2y/p(l-p)(zo + Mp) - (1 - 2p) = 



(13) 



This equation defines the best local basis for the relevant 
Hamiltonian in terms of the same \Gj) and \Ej) states 
defined for the isolated molecule in Eq. © , but with z = 
zo+Mp playing the same role as zq in the definition of p in 
Eq (fTTTf) . With this specific choice of p the uncorrelated gs 
|0) = \G\G2---Gi....) coincides with the gsofthe excitonic 
Hamiltonian, i.e. of H m f + H exc . Since |0) is the gs of 
an Hamiltonian where all on-site terms are retained, it 
also represents the best uncorrelated gs for the system and 
indeed coincides with the mf gs. Within mf the cluster is 
modeled as a collection of non-interacting molecules, each 
molecules experiencing the electric field generated by the 
surrounding^* This local electric field is responsible for 
the renormalization of the energy gap between the local 
I and N states from the value relevant to the isolated 
molecule (2z ) to 2(zq + Mp), relevant to the molecule 
embedded in a cluster of molecules with average polarity 

pM 

The solution of the two-state problem defined by the 
local mf Hamiltonian fully defines the excitonic and ul- 
traexcitonic Hamiltonians. Exciting a molecule costs an 
energy Uojct = 1/ \J p(l ~ p), an d, with p fixed at the mf 
value, H m f can be rewritten as: 



mf 



(14) 



This Hamiltonian counts the number of excited molecules 
in the cluster and assigns energy Hujct to each of them. 
We can also write explicit expressi ons for th e local tran- 
sition dipole moment, per = Po\/ p(l — p), and for the 
local mesomeric dipole moment (the difference between 
the dipole moment in the \Ei) and \Gi)), Ap = p, {l—2p). 
Then, the first term in H ex (Eq. JUJ) describes the inter- 
actions between transition dipole moments on different 
molecules: it corresponds to J-type interactions in the 
standard excitonic model, the only term, beyond H m f in 
Eq. I|14|) , that survives in the excitonic approximation for 
apolar molecules»i* The second term in Eq. (|llfl accounts 
for electrostatic interactions between permanent molec- 
ular dipole moments: it vanishes for apolar molecules. 
H uex in Eq. (|12|) collects ultraexcitonic terms, i.e. terms 
that mix up states with a different number of excitons. 
The first term in Eq. (|12fl has the same origin as the 
exciton hopping, and is usually referred to as the non 
Heitler-London term^a The second term in the ultraex- 
citonic Hamiltonian vanishes for apolar molecules. 

Of course the Hamiltonian in Eq. J§J exactly maps 
on the general Hamiltonian derived by Agranovich for a 
cluster of two-level molecules only interacting via elec- 
trostatic forces^ However, in the standard excitonic ap- 
proximation to the general Hamiltonian, the local ground 
and excited states coincide with the eigenstates of the 
isolated molecule. Environmental contributions to the 
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excitation energies (the term proportional to Mp in the 
first term of H m f) are accounted for, but any mixing 
of the two local eigenstates due to environmental effect 
is neglected (the term proportional to Mp in the sec- 
ond term in H m f is disregarded)^ This amounts to a 
complete neglect of the molecular polarizability: the in- 
teraction with the surrounding only affects excitation en- 
ergies, but not local wavefunctions. Most often the local 
ground and excited states are defined in terms of effec- 
tive molecular states as affected by environmental inter- 
actions, but the parameters entering the excitonic model 
are then freely adjustable and can hardly be related to 
the molecular structure and/or to the supramolecular 
arrangement. Our choice of the local ground and ex- 
cited states as the eigenstates of H m f not only defines 
the best excitonic basis and relates the excitonic and mf 
pictures, but also leads to an unambiguous definition of 
the excitonic and ultraexcitonic Hamiltonians. In fact, 
the mf model relevant to a cluster of non-overlapping 
molecules is defined, for any specific cluster geometry, 
given the model Hamiltonian for the isolated molecule. 8 
At variance with the standard implementation of the ex- 
citonic model, for either apolar or polar molecules, in 
our approach all the parameters entering the excitonic 
(and ultraexcitonic) Hamiltonian, including the excita- 
tion energy, the exciton hopping term and the exciton- 
exciton interaction, are defined by the solution of the 
(self-consistent) mf local Hamiltonian and can be related 
to the molecular and supramolecular structure. The solu- 
tion of the local mf problem is a trivial task for clusters of 
polar and polarizable molecules described in this paper, 
where only two states are enough to capture the (low- 
energy) molecular physics. For MM, like those made up 
of non-polar molecules, where dispersion forces dominate 
intermolecular interactions, setting up and solving the 
local mf problem is a more delicate problem since sev- 
eral molecular states must be explicitely accounted for. 
In any case, uncovering the relation between the local mf 
solution of the cluster Hamiltonian and its excitonic (and 
ultraexcitonic) description is an important result since it 
allows to set up the description of a MM based on infor- 
mation relevant to the isolated molecule, opening the way 
to understand the properties of MM from the molecular 
up to the supramolecular level. 



IV. THE GROUND STATE 

Figure [21 shows the evolution of the gs polarity vs the 
inverse intermolecular distance for A, B, and C clusters 
with v = 1 and two different zq values. In cluster A 
repulsive intermolecular interactions disfavor the charge 
separation on molecular sites and p decreases with w. 
Just the opposite occurs in B and C clusters where at- 
tractive interactions favor charge separation. An isolated 
molecule with an almost neutral gs (p ~ 0.15 at w = for 
Zq = 1) becomes even more neutral in the repulsive, A- 
cluster, but, with increasing intermolecular interactions, 



it becomes more and more I in attractive (B or C) lat- 
tices. For the data in Fig[21 , the N chromophore crosses 
the p = 0.5 interface separating a N from an I gs at w ~ 2 
in C-lattice, i.e. at an inter-chromophore separation ~ 5 
A for typical molecular length ~ 10 A. In B-lattices the 
interface is located at w ~ 0.75, corresponding to inter- 
molecular contacts ~ 3 A. Conversely, molecules with a 
zwitterionic gs in the gas phase (p ~ 0.85 at w = for 
2o=-l) can be driven to a N gs in the repulsive, A-lattice 
at w ~ 1.5, corresponding to intermolecular distances 
~ 7 A. Even more interesting is the observation that 
the same molecule in different supramolecular arrange- 
ments can have qualitatively different gs, and hence dis- 
tinctively different properties. Consider e.g. an isolated 
molecule with an N gs (zo = 1): in an A- type lattice 
with w ~ 3 (i.e. for intermolecular distances ~ one third 
of the molecular length) the molecular polarity stays ba- 
sically unaffected (it slightly decreases down to ~ 0.1). 
The same molecule in a C-lattice with similar intermolec- 
ular distance becomes instead very polar, with p ~ 0.9 
corresponding to an almost zwitterionic gs. Since the 
properties of push-pull chromophores largely depend of 
ptpL very different behavior is predicted for different MM 
made up by the same molecular unit. 

Dashed lines in Fig [21 show mf results for the gs po- 
larity: as far as p is concerned, the mf approach works 
well at least for not too large interactions. Within mf we 
can also understand the qualitatively different behavior 
of attractive and repulsive lattices. As discussed in the 
previous Section, within mf p is given by Eq. J2J), but 
with z = zq + Mp playing the same role of zo, so that: 24 

8zq dz \ dz) ^ ^ 

p(z) has a negative slope, so that for repulsive lattices 
with M > 0, the slope of p(zo) decreases in magni- 
tude with M, justifying the smooth evolution of the p(w) 
curves in FigEK- Just the opposite occurs for attractive 
lattices (M < 0). In this case dp/dzo becomes more 
negative with increasing the strength of inter-site inter- 
actions and, for M < —2, it changes its sign, marking 
the occurrence of a discontinuous phase transition from 
the N to the I regime. The mf N-I crossover is located 
at z = 0, and S-shaped p{w) curves are therefore cal- 
culated within mf for systems with z > 1, as shown in 
Fig. [21 for B and C lattices with v — 2 and zq = 1.5 and 
2. Regions with negative slope in the mf p(w) curves 
correspond to unstable states, and a bistable regime ap- 
pears where two stable states (dp/dw > 0) coexist for 
the same wm^ Demonstrating the appearance of a dis- 
continuous interface from exact diagonalization of finite 
clusters Hamiltonians is hardly possible. In fact only 
stable states are accessible by exact diagonalization: nei- 
ther the unstable nor metastable states appear in exact 
curves. In any case, results from exact diagonalization 
(symbols in Fig. [SJ) show, at large zq, a well pronounced 
discontinuity, consistent with the discontinuous behavior 
recognized within mf. 
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FIG. 2: Ionicity vs the inverse intermolecular distance for 
A, B and C, lattices with v=l and zo = ±1 Dashed lines: mf 
results for N — oo; continuous lines: exact results for TV = 16. 




FIG. 3: Exact ionicity vs the inverse intermolecular distance 
for attractive, B and C, lattices with 7V= 16, v—2 and zo= 
1.5 (squares) and 2 (circles). Dashed line show N = oo mf 
results. 



Isolated (gas phase or solvated) molecules cannot sup- 
port any kind of phase transitions: the appearance of a 
discontinuous N to I crossover is one of the most impres- 
sive evidences of the cooperative behavior of MM made 
up of polar and polarizable chromophores. In the prox- 
imity of the discontinuous crossover large deviations from 
the mf picture are expected on general grounds, due to 
the failure of the non-correlated description of the gs. 
Correlations effects in this region have important spec- 
troscopic consequences and will be discussed in detail in 
Sect. |W| 



V. THE EXCITATION SPECTRUM 

The solution of the local mf problem defines the best 
local basis |Gj) and \Ei) for the excitonic problem and 
hence the optimal excitonic approximation for the Hamil- 
tonian in Eq. Q. Figs. 2IE1 andEl compare the complete 
set of exact eigenstates (circles) with the excitonic eigen- 
states (crosses) for A, B, and C clusters, respectively. To 
avoid overcrowding, data are shown for TV = 6. In the 
figures, the excitation energy (i.e. the energy minus the 
gs energy) is reported for each eigenstate vs n = hi), 
the exciton number. We set zq = —/ + 1 for repul- 
sive/attractive lattices (cf Fig. For each cluster we 
consider two different w values, corresponding to inter- 
actions of weak and medium strength (upper and lower 
panel, respectively). For reference purposes, the inset in 
each panel shows the mf p(w) curve, with the vertical 
dotted line marking the relevant w values. 

Within the excitonic approximation n is conserved: ex- 
citonic eigenstates line up exactly in Figs. [3 andHJ 
defining bands of states with integer n. The excitonic 
bandwidth measures the amount of mixing of states with 
the same n, induced by the exciton hopping (the first 
term in Eq. iJHJ, corresponding to J-like interactions in 
the standard excitonic model) and by the exciton-exciton 
interaction term (the second term in Eq. (|llf> . relevant 
to polar molecules). This last term affects the energy of 
multi-exciton states lowering (increasing) the energy of 
states with several nearby excitations in attractive (re- 
pulsive) lattices. 

The ultraexcitonic Hamiltonian in Eq. (|12f> mixes up 
states with different n. For weak interactions (upper pan- 
els in Fig. 21 [SJ and EJ) the ultraexcitonic mixing is not 
large, and n is approximately conserved also for the exact 
eigenstates (circles). But, for larger interactions (bottom 
panels), the combined role of exciton-exciton interactions 
and of the ultraexcitonic mixing can have important con- 
sequences. 

Figs. 01 and carry one more information: The 
size of the error-bars attached to each symbol are pro- 
portional to the transition dipole moment from the gs 
toward the relevant eigenstate. A general selection rule 
applies to both exact and excitonic eigenstates: only zero 
wavevector (k = 0) states are reached by optical transi- 
tion. Moreover, within the excitonic approximation, only 
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FIG. 4: Excitation spectrum of an A-lattice with N = 6, w=l, 
zq = — 1 and two different w. Circles and crosses show the 
excitation energy for exact and mf eigenstates, respectively, 
against the number of excitons. States on the zero energy 
axis correspond to the gs. Error bars measure the squared 
transition dipole moment from the gs to the relevant states. 
Insets show the p(w) mf curve for the relevant parameters, 
with the dotted vertical line marking the w value for which 
results are reported in the parent panel. 



FIG. 6: The same as in Fig. 2] but for C-lattices with v - 
Zq=1 and two different w. 



1, 




FIG. 5: The same as in Fig. 0] but for B-lattices with v = 1, 
«o=l and two different w. 



one excitation can be created upon one-photon absorp- 
tion so that only the state with k = and n = 1 has a 
finite transition dipole moment from the gs. This state 
is at the top of the n = 1 exciton band in the repulsive 
A lattice (J > 0, H-aggregates) and at the bottom of 
the one-exciton band in the attractive, B-lattice (J < 0, 
J-aggregate). C-lattice with two molecules per unit cell 
is different: the inversion symmetry that relates the two 
molecules in the unit cell allows to classify the two k = 
one-exciton states as symmetric and antisymmetric, and 
only the antisymmetric state is optically active. This 
state corresponds to a k = it state in the extended zone 
representation, and for this attractive lattice (J < 0) the 
optically allowed state actually lies at the top of the one- 
exciton band. 

The ultraexcitonic mixing relaxes the An = 1 selection 
rule of the excitonic model and the oscillator strength can 
spread over several k — states. However, at least for 
the cases shown in Figs. 0] and a single excited 
state retains most of the oscillator strength, leading to a 
spectral behavior qualitatively similar as in the excitonic 
picture. However, for not too weak interactions (bottom 
panels in Figs. 01 [51 and ® the exciton number in the 
most optically allowed state distinctively deviates from 
the excitonic result (n — 1). More generally, the exci- 
tonic description of the photoexcited state is not very 
accurate. Consider, just as an example, a B-lattice with 
parameters as in Fig. \5\ upper panel. As for the data re- 
ported in the figure are concerned, the excitonic picture 
works fairly well. However, the average ionicity of the 
optically allowed state calculated in the excitonic picture 
is about ~ 15% smaller than the exact result. The differ- 
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ence becomes ~ 40% for the parameters in the bottom 
panel of the same figure. The discrepancy between exci- 
tonic and exact descriptions of the states reached upon 
photoexcitation becomes more and more evident as the 
discontinuous charge crossover is approached, and will be 
discussed in greater detail in Sect. IVIII 

VI. STATIC SUSCEPTIBILITIES 

Collective effects that are discussed in the previous Sec- 
tion with reference to excitation spectra, must also show 
up in linear and non-linear susceptibilities of MM. Here 
we focus attention on the linear susceptibility a and on 
the first hyperpolarizability f3, defined by the expansion 
of the gs polarization, P (the dipole moment per unit 
volume) vs a static electric field, F, as follows: 

P(F) = P(0) + aF+ ^(3F 2 + ... (16) 

Continuous lines in Fig. [7| report exact a and (3 val- 
ues calculated as a function of w for 16-site clusters with 
v = 1. Top panels refer to an A-cluster with zq = — 1, 
middle and bottom panels refer to B an C clusters, re- 
spectively, with Zq = 1. The cluster susceptibilities de- 
viate considerably from those relevant to the isolated 
molecules, and a large and non-trivial dependence of the 
optical responses on the supramolecular arrangement is 
observed. This is due to the large response of polarizable 
molecules to the local electric fields in the material, and 
has important, and so far not fully appreciated, conse- 
quences. 

Some degree of non-additivity of the optical responses 
of clusters of push-pull chromophores could have been 
guessed from the large dependence of the gs molecular 
polarity on supramolecular interactions, as shown in Fig. 
121 and |3J For push-pull chromophores in fact linear and 
non-linear optical responses strongly depend on p. Since 
closed expressions for a and f3 as a function of p are 
available^ the simplest approach to static susceptibili- 
ties relies on a mf oriented gas picture. Within mf the 
gs of the cluster describes a collection of non-interacting 
molecules. The responses of a cluster of molecules can 
accordingly be calculated as the sum of contributions 
from the molecules in the relevant mf gs. The suscep- 
tibilities obtained in this approximation (dotted lines in 
Fig. 6) deviates considerably from exact results, a par- 
ticularly impressive result in view of the reliability of mf 
estimates of p for the same clusters (cf Fig. |2J. The 
oriented gas approach to susceptibilities fails since the 
response of a molecule to an applied field strongly de- 
pends on the molecular environment. 

As a matter of fact, the mf approximation leads 
to fairly accurate estimates of the susceptibilities for 
MM, provided the oriented-gas approximation is relaxed. 
Dashed lines in Fig. [7| show the susceptibilities calcu- 
lated as successive P-derivatives of the polarization of 
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FIG. 7: Linear susceptibility (a) and first hyperpolarizabil- 
ity (/?) vs the inverse intermolecular distance for A, B and C 
clusters (upper, middle and lower panels, respectively) with 
TV = 16, v=l and zq — ±1. Dotted lines: mf-oriented 
gas approximation; dot-dashed line: excitonic approximation; 
dashed lines: mf approximation; continuous lines: exact re- 
sults. The j3 response for C-clusters vanishes by symmetry. 



the cluster, calculated in the mf approximation. Apart 
from small deviations that appear just in the crossover 
regime where the mf estimate of p is slightly inaccurate by 
itself (cf Fig. |3J) , the mf susceptibilities nicely agree with 
exact results. The large deviations between the suscepti- 
bilities calculated as sum of molecular mf contributions 
(mf oriented gas approximation, dotted lines in Fig. 0) 
and those obtained from the successive derivatives of the 
mf polarization of the cluster (mf results, dashed lines 
in Fig. UJ are a direct measure of collective effects, and 
can be understood since the polarization of the cluster is 
the sum of the (oriented) dipole moments calculated for 
each molecule, but its derivatives with respect an applied 
field largely deviate from the sum of the corresponding 
derivatives (as imposed in the oriented gas model). 

Sum over state (SOS) expressions directly link static 
susceptibilities to the excitation spectrum^ offering a 
way to estimate the magnitude of excitonic and ultraex- 
citonic contributions to static optical responses. For sure 
susceptibilities evaluated from the successive derivatives 
of the gs polarization and from SOS expressions do co- 
incide, provided they are obtained from the eigenstate of 
the same HamiltonianiSi Then exact SOS susceptibilities 
coincide with the F derivatives of the exact gs polariza- 
tion. Similarly, the susceptibilities obtained in the mf 



oriented gas approach, coincide with the SOS suscepti- 
bilities calculated for a collection of non-interacting (mf) 
molecules. The large deviations of the mf-oriented gas 
susceptibilities from exact results can then be ascribed to 
the combined effect of excitonic and ultraexcitonic mix- 
ing in the excitation spectrum. The complete diagonal- 
ization of the excitonic Hamiltonian (H m f + H ex ) and 
the subsequent calculation of SOS susceptibilities leads 
to results reported as dot-dashed lines in Fig. Ac- 
counting for excitonic interactions considerably improves 
over the oriented gas estimate of static susceptibilities. 
However, sizeable deviations from exact results are ob- 
served, pointing to the importance of ultra-excitonic ef- 
fects. In particular we underline that large ultraexcitonic 
contributions to both a and /3 are observed even in re- 
gions where the mf approximation gives a fairly accurate 
description of both the gs polarity and of the static sus- 
ceptibilities. Neglecting the molecular polarizability in 
the description of excited states, as imposed in the exci- 
tonic approach, is a dangerous approximation for linear 
and non-linear optical responses of molecular materials 
made up of polar and polarizable molecules. 
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VII. COLLECTIVE AND COOPERATIVE 
EFFECTS AT THE DISCONTINUOUS CHARGE 
CROSSOVER 

As discussed in Sect. 11111 the mf approximation ap- 
plied to attractive lattices with z > 1 leads to S-shaped 
p(w) curves (cf FigJSJ, that signal the occurrence of a 
discontinuous phase transition. The excitation spectrum 
of the system near the discontinuous charge crossover is 
interesting and deserves detailed discussion. Figs. IHland 
[5] show the excitation spectrum calculated for B and C 
clusters in the N regime, but very near to the discontin- 
uous crossover. Results are shown for clusters of 10 sites 
and, for the sake of clarity, only zero wavevector eigen- 
states are reported, with circles and crosses referring to 
exact and excitonic results, respectively. 

Near the discontinuous crossover the excitation energy 
(Hujct hi Eq. I|14[l ) is at a minimum and exciton hop- 
ping and exciton-exciton interactions are large enough 
to induce a partial overlap of excitonic bands (cf crosses 
in Figs |H1 13) • Ultraexcitonic mixing is therefore very ef- 
fective: the exciton number is not even approximately 
conserved in exact eigenstates. Moreover, the energy of 
exact excited states show a non-monotonic dependence 
on the exciton number. In particular, in B and C lattices 
the lowest excited state corresponds to a state with more 
than 3 and 2 excitons, respectively. This result contrasts 
sharply with the excitonic picture, where the excitation 
energy increases monotonously with the exciton number, 
so that the lowest excitation energy corresponds to the 
creation of a single exciton. 

As a matter of fact, the excitonic picture is completely 
spoiled in the bistability region. Whereas this is not par- 
ticularly surprising, it is important to realize that many 



FIG. 8: Excitation energies for the zero wavevector eigen- 
states of a 10-site B-cluster with v = 2, zo = 1.5 and w = 0.69 
reported against the number of excitons (upper panels) and 
the number of I molecules (bottom panels). Circles and 
crosses refer to exact and excitonic eigenstates, respectively. 
States on the zero energy axis correspond to the gs. Error 
bars measure the squared transition dipole moment from the 
gs to the relevant states. The inset shows the correspond- 
ing p(w) mf curve, with the vertical dotted line marking the 
relevant w value. 



concepts usually adopted to discuss excitations in molec- 
ular crystals fail there. In particular, the very same con- 
cept of excitation is challenged by exact results for B and 
C clusters near the discontinuous crossover. Both lattices 
in fact have in the gs a number of excitons distinctively 
larger than zero (cf Figs|Hland[!||where the exact gs eigen- 
state has n > 0) . In these conditions the definition of a 
local excitation is somewhat artificial. We therefore go 
back to the original \DA) and \D + A~) basis and, in the 
lower panels in Figs|S| andEjwe report the same excita- 
tion spectrum as in the upper panels, but with the ab- 
scissa axis now measuring the average number of fully-I 
(\D + A~)) sites, i.e. the expectation value of hi = J2i Pi- 
The operator hi is diagonal on the excitonic eigenstates: 
in the excitonic approximation nj and n are related by a 
simple expression, ni = Np + (1 — 2p)n, where ni, and 
n are the expectation values of the corresponding opera- 
tors in the given state and p is the mf molecular ionicity. 
Excitonic eigenstates then line up vertically also when 
reported against ni (cf crosses in Figs. QOandEJ). For 
exact eigenstates instead n and ni carry independent in- 
formation. 
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FIG. 9: The same as in Fig. |H|for a C cluster with v = 2, 
20 = 1.5 and w = 1.7. For graphical purposes the error bars 
measuring the transition dipole moments have been magnified 
by a factor of two with respect to those reported in Fig. |H] 



In any case, ni is finite in both the exact and exci- 
tonic gs to indicate a finite p. The difference between 
the excitonic and exact ni of course parallels the size- 
able difference between the mf and exact p observed near 
the discontinuous crossover (cf Fig. |3J). More interesting 
is the behavior of the lowest excited state: in both B and 
C lattices the lowest excitation corresponds to a state 
whose ni appreciably differs from the gs value: the low- 
est exact excited states has a very different nature from 
the gs. This contrasts with the excitonic picture where 
the lowest excited state is a state with a single excitation, 
whose ionicity differs from the gs ionicity by |1 — 2p\ <1. 

Attractive supramolecular interactions in B and C lat- 
tices lead to a qualitatively similar p{w) dependence (cf 
Figs [21 EJ). The two lattices however have different sym- 
metry and hence a qualitatively different spectroscopic 
behavior, as already discussed in Sec. E] Lattice B has 
one molecule per unit cell and all k = state are optically 
allowed, even if most of them have negligible intensity. 
The two molecules per cell in lattice C are exchanged by 
reflection, so that only k = antisymmetric states are 
accessible by one-photon absorption from the (totally- 
symmetric) gs: about half of the exact k = states in 
Fig. EJare actually dark states. The special low-energy 
state with large n is a dark state in C-clusters, but corre- 
sponds to the state with the largest oscillator strength in 
B-clusters (cf Fig. 8, and 9). The spectroscopic behavior 
of B-lattices is therefore particularly interesting and will 



be discussed in greater detail below. 

Data in Fig. [8] show that upon absorption of a sin- 
gle photon about 3 excitons are created in B-lattice, or 
~ 4 molecules are switched from N to I. This implies of 
course that the motion of ~ 4 electrons is driven by the 
absorption of a single photon. In B-lattices near the dis- 
continuous interface the primary photoexcitation event 
therefore corresponds to a multielectron transfer^ This 
sharply contrasts with the excitonic picture where a sin- 
gle molecule is excited by a single photon, leading to the 
transfer of no more than a single electron (specifically 
1 1 — 2yo| < 1 electrons are transferred upon photoexcita- 
tion). Multielectron transfer is a new phenomenon with 
no counterpart in the standard description of optical exci- 
tations in MM. The theoretical and practical implications 
of multielectron transfer arc hardly overemphasized. 29 
Here we underline that so far multielectron transfer was 
discussed as a secondary photoexcitation event: the ab- 
sorption of a photon induces a single electron transfer, 
then cooperative interaction with slow degrees of freedom 
(either vibrational or environmental) leads to multiple 
electron transfer via a cascading effecti^S In our model 
direct (vertical) photoexcitation of several electrons is 
made possible by the correlation of electronic motion on 
different (non-overlapping) molecules for systems at the 
discontinuous charge crossover^ 

Data in Fig. [S] refer to a 10 site cluster, but multi- 
electron transfer safely survives in longer chains. With 
reference to the excited state with the largest transition 
dipole moment, we calculate both the excitation energy 
and the average number of molecular sites turned I upon 
excitation (vi, defined as the difference between ni in the 
excited and ground state). Left panels in Fig. 1101 show 
the N-dependence of these two quantities calculated for 
a B-lattice with v — 2, zq = 1.5, and several w. For w 
increasing from 0.60 to 0.69 the system is driven towards 
the interface (cf Fig. |3J) and vj smoothly extrapolates to 
N — > oo, with values that increase from ~ 1.2 (w — 0.60) 
to 2.5 (w—0.68). Finite-size effects are larger at w=0.69, 
but vi clearly extrapolates to even larger (~3) values 
there. The corresponding excitation energy smoothly 
extrapolates towards finite values in the N — > oo limit 
(Fig. 1101 upper panel): in an infinite linear chain with 
B-geometry, near the discontinuous charge crossover the 
absorption of a single optical photon (with energies ~ 1 
eV for the parameters in Fig. 1101 and \/2t=leV) drives 
the coherent motion of several (~ 3) electrons. 

The right panels in Fig. EH show the behavior of a 
B-lattice near to a discontinuous charge crossover with 
a much larger ionicity jump with respect to case shown 
in the left panel (cf. Fig. |3J). Results are even more 
impressive there: depending on the model parameters, in 
fact up to 8 electrons can be transferred by absorption of 
a single photon. However, the finite-size analysis requires 
in this case much longer chains than addressed here (at 
least a few times longer than the number of electrons 
transferred). 

To better understand the physics of multielectron 
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FIG. 10: The 1/iV dependence of the excitation energy and 
of the average number of molecules turned I upon photoexci- 
tation. Results refer to the state with the largest transition 
dipole moment for B-clusters with v — 2, and zo — 1.5, and 2 
(left and right panels, respectively). The numbers labelling 
each curve show the corresponding w value. Lines are drawn 
as guide for eyes. 



transfer we define the following Z-th order correlation 
function: 



/v 



N 



iPi+l---Pi+l-l) 



(17) 



i=l 



fi vanishes exactly for uncorrelated states, that is for 
states that can be reduced to the product of local molec- 
ular states. Positive (negative) /; indicates an increased 
(decreased) probability of finding I nearby I molecules 
with respect to the uncorrelated state with the same av- 
erage polarity. The upper panel of Fig. II II shows the l- 
dependence of // calculated for the most optically active 
state of a B-cluster with the same parameters as in Fig. 
1101 just at the N-I crossover (w — 0.69). In this state 
a non-negligible weight is found of wavefunctions with 
several (2-4) nearby fully-I molecules. Then the photo- 
induced multielectron transfer corresponds to a concerted 
electronic motion occurring on several nearby molecules. 

The bottom panel in Fig. ^] shows the correlation 
function calculated for the gs: quite predictably the over- 
all degree of correlation is smaller in the ground than 
in the excited state (notice the different scale on the y- 
axis in the two panels in Fig. Ill|) . In any case data in 
Fig. ^] demonstrate a finite contribution to the gs of 
I-droplet states, i.e. of states with two or more nearby 
I-sites. Near the discontinuous charge crossover the mf 
description of the gs fails: the exact gs is not the vacuum 
state (cf the finite n e relevant to the exact gs in Fig. [SJ 
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FIG. 11: The / dependence of /; correlation function calcu- 
lated for the ground state and for the excited state with the 
largest transition dipole moment for a B-lattice with v = 2, 
zo = 1.5 and w = 0.69. Triangles, squares and circles shows 
results for iV= 12, 14, and 16, respectively. 



and cannot be described as the direct product of local 
molecular states (cf the finite /; in the lowest panel in 
fig. Ill|) . Within the proposed model for non-overlapping 
molecules classical electrostatic interactions lead, near 
the discontinuous crossover, to an intrinsically collective 
and correlated zero-temperature gs. The finite amplitude 
of I-droplet states in this gs is the key to understand mul- 
tielectron transfer: in B-clusters, I droplet states have a 
finite (and large) permanent dipole moment: their finite 
amplitude in the gs is therefore responsible for the ap- 
pearance of sizeable transition dipole moments from the 
gs towards states with a large I-droplet character. This 
observation also sheds light on the qualitatively different 
behavior of C and B lattices. For sure the discontinu- 
ous crossover in C-lattices is similar in nature to that 
observed in B-lattices and correlated ground and excited 
states are expected for these lattices too. However C- 
lattices are centrosymmetric: droplet states do not have 
permanent dipole moments and do not contribute to the 
transition dipole moments. 

The behavior of systems lying near the discontinuous 
crossover, but on the I side is similar to that of N systems, 
described above, provided the role of \DA) and \D + A~) 
is interchanged. Similar behavior as described in Figs. 
1101 and ^2 can be obtained in fact if uj, the number 
of sites turned I upon photoexcitation is substituted by 
vjsr, the number of sites turned N. Similarly 1 — pi must 
substitute pi to define a correlation function measuring 
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the probability, in each relevant state, of droplets of I 
nearby N-molcculcs. 



VIII. DISCUSSION 

In this paper we have presented a simple and interest- 
ing model for clusters of polar and polarizable molecules. 
Each molecular site is described in terms of a two-state 
model. Only classical electrostatic intermolecular inter- 
actions are accounted for, neglecting any intermolecu- 
lar overlap of wavefunctions on different molecules. The 
model applies quite naturally to MM (crystals, aggre- 
gates, films..) based on push-pull chromophores, an in- 
teresting class of molecules quite extensively investigated 
for NLO applications. 

Push-pull chromophores are largely polar and polariz- 
able molecules, and electrostatic interactions can only be 
understood within models that properly account for the 
molecular polarizability at all orders: each molecule in 
fact experiences in the material the local electric fields 
generated by the surrounding molecules in a non-trivial 
feed-back mechanism that is responsible for large col- 
lective and cooperative effects. The Mulliken model, 
adopted for the molecular sites, is particularly interest- 
ing in this respect since it accounts for the molecular 
polarizability and hyperpolarizabilities, 25 in a two-state 
approach that represents the simplest picture to describe 
polar molecules. 

Our model of course maps exactly onto the Hamil- 
tonian derived several years ago by Agranovich for a 
molecular crystal described as a collection of two-state 
molecules interacting via electrostatic forces^ Being in- 
terested in the very weak interaction limit, Agranovich 
defined the excitonic Hamiltonian by disregarding all 
terms in the complete Hamiltonian not conserving the 
exciton number i— For polarizable molecules we propose 
a different strategy: via a rotation of the local (on-site) 
basis we retain the local term not conserving n. In other 
words we recognize that the optimal local basis for the 
excitonic problem coincides with the eigenstates of the lo- 
cal self-consistent mf problem. The relation between the 
mf and the excitonic approach is fairly fundamental. In 
the excitonic model excitations are created on top of an 
uncorrelated gs: the mf solution then gives the best gs for 
the excitonic problem. The definition of the local states 
for the excitonic problem as the eigenstates of the lo- 
cal mf Hamiltonian unambiguously defines the excitonic 
and ultraexcitonic Hamiltonian for any supramolecular 
arrangement, based on the adopted model for the iso- 
lated molecule. Uncovering a direct link from the molec- 
ular to the supramolecular description gives an important 
contribution to our understanding of MM, and is a funda- 
mental step to devise approaches to guide the chemical 
synthesis of functional MM from the molecular to the 
supramolecular level. 

The mf description of MM can be set up for dif- 
ferent kinds of molecules described at different level 



of sophistication. 8 So our approach is widely applica- 
ble. Here we have discussed its application to a simple 
model for clusters of polar and polarizable chromophores 
where particularly important collective and cooperative 
effects are expected. Collective effects show up in the gs 
with the large dependence of the molecular polarity on 
supramolecular interactions. The effect of the surround- 
ing medium, and specifically of the solvent, in tuning 
the polarity of push-pull chromophores has already been 
underlinedj22i2i and is a natural consequence of the large 
molecular polarizability. Here we extend this concept 
to supramolecular interactions in MM and demonstrate 
that they can widely tune the molecular polarity, and 
eventually drive it across the N-I interface. The major 
difference with respect to solvated molecules is the pos- 
sibility for MM to support true phase transitions, as the 
extreme manifestation of cooperative behavior. 

The properties of isolated push-pull chromophores, and 
most interestingly, their NLO responses strongly depend 
on the molecular polarity^ the large variation of p in- 
duced by supramolecular interactions then immediately 
suggests that NLO responses are strongly affected by 
supramolecular interactions. Apart from a narrow re- 
gion around the crossover regime, the mf approach nicely 
approximate the NLO responses of the material. Huge 
collective effects are demonstrated that are only partly 
accounted for within the excitonic approximation. The 
large ultraexcitonic corrections in Fig. Qsuggest that the 
excitonic model does not provide a reliable approxima- 
tion scheme for the calculation of NLO responses of MM 
based on largely polarizable molecules. 

The traditional strategy to get molecular materials 
with large (3 responses relies on the optimization of /3 
at the molecular level, in the tacit assumption that ma- 
terials based on molecules with the largest (3 will have the 
largest responses. But, in agreement with experimental 
observation;22i2&2i results in Fig. [7| demonstrate that 
the properties of MM based on push-pull chromophores 
are far from additive. To guide the synthesis of opti- 
mized materials, structure properties relationships must 
be devised and fully understood not only at the molecular 
level but also at the supramolecular level. The detailed 
analysis of the simple model for supramolecular clusters 
presented here is a first step in this direction. 

Deviations of the excitation spectrum from the ex- 
citonic description are predictably sizeable for materi- 
als with medium-large supramolecular interactions, with 
more important effects for materials located near the N-I 
interface where the molecular polarizability is at a max- 
imum. More interesting is the observation of new phe- 
nomena driven by cooperativity near the discontinuous 
charge-crossover of attractive lattices. The lowest excita- 
tion of a N (I) lattice near the N-I discontinuous crossover 
creates an droplet of I (N) molecules whose size increases 
as the system is driven towards more discontinuous inter- 
faces. These droplet states, that appear as low-lying ex- 
citations near the discontinuous N-I crossover, share the 
same physics with the charge-density-wave droplets ob- 
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served in half-filled extended Hubbard ID chains at the 
discontinuous phase transition from the Mott-insulator 
to the charge-density wave stated and are due to the 
nucleation of I (N) domains in the presence of large in- 
tersite interactions. The main novelty of our result in 
this respect is the observation that these droplet states 
are, in a non-centrosymmetric environment, optically al- 
lowed, so that the absorption of a single photon in a N 
(I) material can directly create an I (N) droplet extend- 
ing over several molecules. The primary photoexcitation 
event is in these conditions a multielectron transfer that 
is a direct consequence of the correlation of the electronic 
motion among different non-overlapping molecules. 

We have discussed the gs and the excitation spec- 
trum of a model for polar and polarizable molecular 
units, only interacting via electrostatic forces. The 
model naturally applies to MM based on push-pull chro- 
mophores, but also offers a rough description of organic 
CT salts with a mixed stack motif. In these materials, 
7r-conjugated molecules with a strong electron D and A 
character alternate in a one-dimensional stack: the size- 
able overlap between adjacent molecules leads to charge- 
resonance and hence to fractional charges on molecu- 
lar sites. The N-I phase transition observed in these 
materials^ is a fascinating phenomenon whose physics 
is actively investigated i^l In the early days of the N-I 
transition, it was modeled by reducing the stack of over- 
lapping molecules into a collection of non-overlapping 
DA pairSfSi i.e. by adopting exactly the same model 
proposed here for B-lattice. Of course this approxima- 
tion is very rough and more refined models are needed 
to understand the complex and variegated behavior of 



CT salts with a mixed stack motif However, the N to 
I crossover discussed here for attractive lattices of po- 
lar and polarizable molecules shares the same physics as 
the N-I transition in CT salts. In this respect we un- 
derline that anomalous spectral features for CT salts in 
the region of the discontinuous N-I transition were al- 
ready discussed in Ref.391 In particular the importance 
of I/N droplet states as well as the failure of the exci- 
tonic approximation in the proximity of the discontin- 
uous crossover were addressed. Multielectron transfer, 
described here for molecular clusters, indeed offers sug- 
gestive hints on the physics of photoinduced N-I phase 
transitions in CT salts. 28 However CT salts are extended 
systems, with electrons delocalized all along the stack. 
On the opposite, the molecular clusters described here 
are fully localized systems: electrons cannot escape from 
their parent molecule. In this respect we underline that 
the main novelty and interest of the present work is in 
the observation of cooperative phenomena and correla- 
tion effects among non-overlapping molecular units. 
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